Phylogenomics and the rise of the angiosperms.

Angiosperms are the cornerstone of most terrestrial ecosystems and human livelihoods1,2. A robust understanding of angiosperm evolution is required to explain their rise to ecological dominance. So far, the angiosperm tree of life has been determined primarily by means of analyses of the plastid genome3,4. Many studies have drawn on this foundational work, such as classification and first insights into angiosperm diversification since their Mesozoic origins5-7. However, the limited and biased sampling of both taxa and genomes undermines confidence in the tree and its implications. Here, we build the tree of life for almost 8,000 (about 60%) angiosperm genera using a standardized set of 353 nuclear genes8. This 15-fold increase in genus-level sampling relative to comparable nuclear studies9 provides a critical test of earlier results and brings notable change to key groups, especially in rosids, while substantiating many previously predicted relationships. Scaling this tree to time using 200 fossils, we discovered that early angiosperm evolution was characterized by high gene tree conflict and explosive diversification, giving rise to more than 80% of extant angiosperm orders. Steady diversification ensued through the remaining Mesozoic Era until rates resurged in the Cenozoic Era, concurrent with decreasing global temperatures and tightly linked with gene tree conflict. Taken together, our extensive sampling combined with advanced phylogenomic methods shows the deep history and full complexity in the evolution of a megadiverse clade.

Light-sensitive short hypocotyl genes confer symbiotic nodule identity in the legume Medicago truncatula.

Legumes produce specialized root nodules that are distinct from lateral roots in morphology and function, with nodules intracellularly hosting nitrogen-fixing bacteria. We have previously shown that a lateral root program underpins nodule initiation, but there must be additional developmental regulators that confer nodule identity. Here, we show two members of the LIGHT-SENSITIVE SHORT HYPOCOTYL (LSH) transcription factor family, predominantly known to define shoot meristem complexity and organ boundaries, function as regulators of nodule organ identity. In parallel to the root initiation program, LSH1/LSH2 recruit a program into the root cortex that mediates the divergence into nodules, in particular with cell divisions in the mid-cortex. This includes regulation of auxin and cytokinin, promotion of NODULE ROOT1/2 and Nuclear Factor YA1, and suppression of the lateral root program. A principal outcome of LSH1/LSH2 function is the production of cells able to accommodate nitrogen-fixing bacteria, a key feature unique to nodules.

NODULE INCEPTION Recruits the Lateral Root Developmental Program for Symbiotic Nodule Organogenesis in Medicago truncatula.

To overcome nitrogen deficiencies in the soil, legumes enter symbioses with rhizobial bacteria that convert atmospheric nitrogen into ammonium. Rhizobia are accommodated as endosymbionts within lateral root organs called nodules that initiate from the inner layers of Medicago truncatula roots in response to rhizobial perception. In contrast, lateral roots emerge from predefined founder cells as an adaptive response to environmental stimuli, including water and nutrient availability. CYTOKININ RESPONSE 1 (CRE1)-mediated signaling in the pericycle and in the cortex is necessary and sufficient for nodulation, whereas cytokinin is antagonistic to lateral root development, with cre1 showing increased lateral root emergence and decreased nodulation. To better understand the relatedness between nodule and lateral root development, we undertook a comparative analysis of these two root developmental programs. Here, we demonstrate that despite differential induction, lateral roots and nodules share overlapping developmental programs, with mutants in LOB-DOMAIN PROTEIN 16 (LBD16) showing equivalent defects in nodule and lateral root initiation. The cytokinin-inducible transcription factor NODULE INCEPTION (NIN) allows induction of this program during nodulation through activation of LBD16 that promotes auxin biosynthesis via transcriptional induction of STYLISH (STY) and YUCCAs (YUC). We conclude that cytokinin facilitates local auxin accumulation through NIN promotion of LBD16, which activates a nodule developmental program overlapping with that induced during lateral root initiation.

Dominant integration locus drives continuous diversification of plant immune receptors with exogenous domain fusions.

BACKGROUND: The plant immune system is innate and encoded in the germline. Using it efficiently, plants are capable of recognizing a diverse range of rapidly evolving pathogens. A recently described phenomenon shows that plant immune receptors are able to recognize pathogen effectors through the acquisition of exogenous protein domains from other plant genes. RESULTS: We show that plant immune receptors with integrated domains are distributed unevenly across their phylogeny in grasses. Using phylogenetic analysis, we uncover a major integration clade, whose members underwent repeated independent integration events producing diverse fusions. This clade is ancestral in grasses with members often found on syntenic chromosomes. Analyses of these fusion events reveals that homologous receptors can be fused to diverse domains. Furthermore, we discover a 43 amino acid long motif associated with this dominant integration clade which is located immediately upstream of the fusion site. Sequence analysis reveals that DNA transposition and/or ectopic recombination are the most likely mechanisms of formation for nucleotide binding leucine rich repeat proteins with integrated domains. CONCLUSIONS: The identification of this subclass of plant immune receptors that is naturally adapted to new domain integration will inform biotechnological approaches for generating synthetic receptors with novel pathogen "baits."

Genome dynamics explain the evolution of flowering time CCT domain gene families in the Poaceae.

Numerous CCT domain genes are known to control flowering in plants. They belong to the CONSTANS-like (COL) and PREUDORESPONSE REGULATOR (PRR) gene families, which in addition to a CCT domain possess B-box or response-regulator domains, respectively. Ghd7 is the most recently identified COL gene to have a proven role in the control of flowering time in the Poaceae. However, as it lacks B-box domains, its inclusion within the COL gene family, technically, is incorrect. Here, we show Ghd7 belongs to a larger family of previously uncharacterized Poaceae genes which possess just a single CCT domain, termed here CCT MOTIF FAMILY (CMF) genes. We molecularly describe the CMF (and related COL and PRR) gene families in four sequenced Poaceae species, as well as in the draft genome assembly of barley (Hordeum vulgare). Genetic mapping of the ten barley CMF genes identified, as well as twelve previously unmapped HvCOL and HvPRR genes, finds the majority map to colinear positions relative to their Poaceae orthologues. Combined inter-/intra-species comparative and phylogenetic analysis of CMF, COL and PRR gene families indicates they evolved prior to the monocot/dicot divergence ∼200 mya, with Poaceae CMF evolution described as the interplay between whole genome duplication in the ancestral cereal, and subsequent clade-specific mutation, deletion and duplication events. Given the proven role of CMF genes in the modulation of cereals flowering, the molecular, phylogenetic and comparative analysis of the Poaceae CMF, COL and PRR gene families presented here provides the foundation from which functional investigation can be undertaken.

The molecular basis for venation patterning of pigmentation and its effect on pollinator attraction in flowers of Antirrhinum.

Pigment stripes associated with veins (venation) is a common flower colour pattern. The molecular genetics and function of venation were investigated in the genus Antirrhinum, in which venation is determined by Venosa (encoding an R2R3MYB transcription factor). Pollinator preferences were measured by field tests with Antirrhinum majus. Venosa function was examined using in situ hybridization and transient overexpression. The origin of the venation trait was examined by molecular phylogenetics. Venation and full-red flower colouration provide a comparable level of advantage for pollinator attraction relative to palely pigmented or white lines. Ectopic expression of Venosa confers pigmentation outside the veins. Venosa transcript is produced only in small areas of the corolla between the veins and the adaxial epidermis. Phylogenetic analyses suggest that venation patterning is an ancestral trait in Antirrhinum. Different accessions of three species with full-red pigmentation with or without venation patterning have been found. Epidermal-specific venation is defined through overlapping expression domains of the MYB (myoblastoma) and bHLH (basic Helix-Loop-Helix) co-regulators of anthocyanin biosynthesis, with the bHLH providing epidermal specificity and Venosa vein specificity. Venation may be the ancestral trait, with full-red pigmentation a derived, polyphyletic trait. Venation patterning is probably not fixed once species evolve full-red floral pigmentation.

Comparative genomics of flowering time pathways using Brachypodium distachyon as a model for the temperate grasses.

Brachypodium distachyon (Brachypodium) is a model for the temperate grasses which include important cereals such as barley, wheat and oats. Comparison of the Brachypodium genome (accession Bd21) with those of the model dicot Arabidopsis thaliana and the tropical cereal rice (Oryza sativa) provides an opportunity to compare and contrast genetic pathways controlling important traits. We analysed the homologies of genes controlling the induction of flowering using pathways curated in Arabidopsis Reactome as a starting point. Pathways include those detecting and responding to the environmental cues of day length (photoperiod) and extended periods of low temperature (vernalization). Variation in these responses has been selected during cereal domestication, providing an interesting comparison with the wild genome of Brachypodium. Brachypodium Bd21 has well conserved homologues of circadian clock, photoperiod pathway and autonomous pathway genes defined in Arabidopsis and homologues of vernalization pathway genes defined in cereals with the exception of VRN2 which was absent. Bd21 also lacked a member of the CO family (CO3). In both cases flanking genes were conserved showing that these genes are deleted in at least this accession. Segmental duplication explains the presence of two CO-like genes in temperate cereals, of which one (Hd1) is retained in rice, and explains many differences in gene family structure between grasses and Arabidopsis. The conserved fine structure of duplications shows that they largely evolved to their present structure before the divergence of the rice and Brachypodium. Of four flowering-time genes found in rice but absent in Arabidopsis, two were found in Bd21 (Id1, OsMADS51) and two were absent (Ghd7, Ehd1). Overall, results suggest that an ancient core photoperiod pathway promoting flowering via the induction of FT has been modified by the recruitment of additional lineage specific pathways that promote or repress FT expression.

IT3F: a web-based tool for functional analysis of transcription factors in plants.

A web-based tool, the Interspecies Transcription Factor Function Finder (IT3F), has been developed to display both evolutionary gene relationships and expression data for plant transcription factors, focussing primarily on the R2R3MYB gene subfamily for proof of concept. The graphical display of information allows users to make direct comparisons between structurally related genes and to identify those genes that are potentially orthologous, thereby assisting with their understanding of gene function. A key feature of the website is the provision of an interrogative phylogenetic tree that allows submission of new sequences corresponding to a transcription factor family or subfamily and maps their relative positions to the products of other genes on an 'existing' tree containing proteins encoded by Arabidopsis and rice genes, along with key proteins encoded by genes from other species that have been characterised functionally. In addition, a feature to select clusters of related sequences has been developed so that more detailed phylogenetic analysis can be performed to highlight potential orthologous and paralogous genes within related clusters. Arabidopsis genes that reside on duplicated regions of the genome are indicated on the tree, providing further information for interpreting gene function. An additional feature of the website allows a selected number of key Arabidopsis and rice microarray experiments to be visualised alongside the tree as a tabulated heat map of expression intensity values. Through this display, it is possible to observe relative expression levels across a whole gene family and the extent to which the expression of closely related genes within subgroups has altered since their ancestral divergence. The website is available at http://jicbio.nbi.ac.uk/IT3F/.

Environmental deterioration increases tadpole vulnerability to predation.

Human-induced environmental change is occurring at an unprecedented rate and scale. Many freshwater habitats, in particular, have been degraded as a result of increased salinity. Little is known about the effects of anthropogenic salinization on freshwater organisms, especially at sublethal concentrations, where subtle behavioural changes can have potentially drastic fitness consequences. Using a species of Australian frog (Litoria ewingii), we experimentally examined the effects of salinization on tadpole behaviour and their vulnerability to a predatory dragonfly nymph (Hemianax papuensis). We found that tadpoles exposed to an ecologically relevant concentration of salt (15% seawater, SW) were less active than those in our freshwater control (0.4% SW). Tadpoles in elevated salinity also experienced a higher risk of predation, even though the strike rate of the predator did not differ between salt and freshwater treatments. In a separate experiment testing the burst-speed performance of tadpoles, we found that tadpoles in saltwater were slower than those in freshwater. Thus, it would appear that salt compromised the anti-predator response of tadpoles and made them more susceptible to being captured. Our results demonstrate that environmentally relevant concentrations of aquatic contaminants can, even at sublethal levels, severely undermine the fitness of exposed organisms.

The basic helix-loop-helix transcription factor family in plants: a genome-wide study of protein structure and functional diversity.

Basic helix-loop-helix (bHLH) transcription factors (TFs) belong to a family of transcriptional regulators present in three eukaryotic kingdoms. Many different functions have been identified for these proteins in animals, including the control of cell proliferation and development of specific cell lineages. Their mechanism for controlling gene transcription often involves homodimerization or heterodimerization. In plants, little is known about the bHLH family, but we have determined that there are 133 bHLH genes in Arabidopsis thaliana and have confirmed that at least 113 of them are expressed. The AtbHLH genes constitute one of the largest families of transcription factors in A. thaliana with significantly more members than are found in most animal species and about an equivalent number to those in vertebrates. Comparisons with animal sequences suggest that the majority of plant bHLH genes have evolved from the ancestral group B class of bHLH genes. By studying the AtbHLH genes collectively, twelve subfamilies have been identified. Within each of these main groups, there are conserved amino acid sequence motifs outside the DNA binding domain. Potential gene redundancy among members of smaller subgroups has been analyzed, and the resulting information is presented to provide a simplified visual interpretation of the gene family, identifying related genes that are likely to share similar functions. Based on the current characterization of a limited number of plant bHLH proteins, we predict that this family of TFs has a range of different roles in plant cell and tissue development as well as plant metabolism.

Transcripts of Vp-1 homeologues are misspliced in modern wheat and ancestral species.

The maize (Zea mays) Viviparous 1 (Vp1) transcription factor has been shown previously to be a major regulator of seed development, simultaneously activating embryo maturation and repressing germination. Hexaploid bread wheat (Triticum aestivum) caryopses are characterized by relatively weak embryo dormancy and are susceptible to preharvest sprouting (PHS), a phenomenon that is phenotypically similar to the maize vp1 mutation. Analysis of Vp-1 transcript structure in wheat embryos during grain development showed that each homeologue produces cytoplasmic mRNAs of different sizes. The majority of transcripts are spliced incorrectly, contain insertions of intron sequences or deletions of coding region, and do not have the capacity to encode full-length proteins. Several VP-1-related lower molecular weight protein species were present in wheat embryo nuclei. Embryos of a closely related tetraploid species (Triticum turgidum) and ancestral diploids also contained misspliced Vp-1 transcripts that were structurally similar or identical to those found in modern hexaploid wheat, which suggests that compromised structure and expression of Vp-1 transcripts in modern wheat are inherited from ancestral species. Developing embryos from transgenic wheat grains expressing the Avena fatua Vp1 gene showed enhanced responsiveness to applied abscisic acid compared with the control. In addition, ripening ears of transgenic plants were less susceptible to PHS. Our results suggest that missplicing of wheat Vp-1 genes contributes to susceptibility to PHS in modern hexaploid wheat varieties and identifies a possible route to increase resistance to this environmentally triggered disorder.

The feeding behaviour of a sit-and wait-predator, Ranatra dispar (Heteroptera: Nepidae): optimal foraging and feeding dynamics.

The feeding behaviour of R. dispar was examined with respect to the proportion of prey contents used, the time between successive captures (intercatch interval) and the feeding time. The feeding process consisted of three stages. (1) Injection of venom, (2) breakdown of tissue/digestive stage and (3) extraction of food. The rate of extraction from an individual prey decreases as its contents are depleted, but was shown to increase significantly during the first 15 min before decreasing. Even after 30 min the extraction rate was still marginally higher than the initial extraction rate. This phenomenon is quite different to what has previously been reported for sucking bugs.There was a negative relationship between increasing prey density and prey depletion, with the predators being significantly more 'wasteful', (i.e. prey were discarded before all extractable food was removed) at the two higher prey densities compared with those at the lower densities. As the prey density decreased from 60 to 1 prey per container, so the resultant intercatch interval, feeding time, and the average dry weight extracted per prey increased. No correlation was found between individual intercatch interval and subsequent feeding time when examined throughout a sequence of eight captures. This is taken to support the optimal feeding model in which the predator reacts to the average profitability of the environment (i.e. mean intercatch interval) rather than as reflected by the amount of food in the gut. The effect of the changing rate of extraction of food during a meal allows Ranatra, when exposed to high prey density, to feed for less than half the time on each prey item that it spends at low densities, and yet still obtain 60% of available food.